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Abstract. The voter model is a paradigm of ordering dynamics. At each time step, 
a random node is selected and copies the state of one of its neighbors. Traditionally, 
this state has been considered as a binary variable. Here, we relax this assumption 
and address the case in which the number of states is a parameter that can assume 
any value, from 2 to oo, in the thermodynamic limit. We derive mean- field analytical 
expressions for the exit probability and the consensus time for the case of an arbitrary 
number of states. We then perform a numerical study of the model in low dimensional 
lattices, comparing the case of multiple states with the usual binary voter model. Our 
work generalizes the well-known results for the voter model, and sheds light on the 
role of the so far almost neglected parameter accounting for the number of states. 



PACS numbers: 89.65.-s, 05.40.-a, 89.75.-k 



Ordering dynamics of the multi-state voter model 



2 



1. Introduction 

Models of ordering dynamics have since long been considered as paradigms of opinion 
dynamics and consensus formation in social systems (JTj). Most of them share the 
fundamental feature that order results from the self-organization of local and usually 
short-range pairwise interactions between agents, a feature that is well exemplified in 
the simplest and most analyzed of them, the so-called voter model (J2j [3j) . In its basic 
formulation, the voter model is defined as follows: Each individual in a population 
(agent) is endowed with a binary spin variable, representing two opposed opinions, and 
taking values a — ±1. At each time step, an agent i, together with one nearest neighbor 
j, are selected at random, and the state of the system is updated as cxj := <Tj, the first 
agent copying the opinion of its neighbor. Starting from a disordered initial state, this 
dynamics leads in finite systems to a uniform state with all individuals sharing the same 
opinion (the so-called consensus). 

From the point of view of social dynamics, the interest in this kind of models is 
mainly focused on the way in which consensus is reached. The approach to this state 
is characterized in terms of the exit probability E(x) and the consensus time T^(x), 
defined as the probability that the final state corresponds to all agents in the state +1 
and the average time needed to reach consensus in a system of size N, respectively, 
when starting from a homogeneous initial condition with a fraction x of agents in state 
+ 1 (1). Due to its simplicity, the voter model dynamics can be exactly solved in regular 
lattices for any number of dimensions (jUEJ). Thus, considering the average conservation 
of magnetization m = X^i "*/^ ^ can De shown that the exit probability is always 
a linear function, E(x) = x. On the other hand, the consensus time starting from 
the homogeneous symmetric initial condition x = 1/2 scales with system size iV as 
T N (l/2) ~ N cS , with N eS ~ N 2 in d = 1, N eS ~ JVlogiV in d = 2, and N eS ~ iV in 
d > 2 (at the mean-field level) (0). Finally, the dependence of consensus time with the 
initial density of +1 spins, starting from homogeneous initial conditions, takes the form 



for d > 2 (|SJ. 

Different variants of the voter model have been considered in the past, including 
the presence of quenched disorder in the form of "zealots" which do not change opinion 
<m IE]), memory and noise reduction ([9]), inertia (flOl) . non-conservative voters ( 1TT1) . non- 
linear interactions (fT2llT3j) . etc.; see Ref. (jl]) for an extended bibliography on this subject. 
A variant that has been considered in several contexts is the multi-state voter model, 
in which each agent can be in one of S different exclusive states or opinions, in analogy 
of the Potts model (I14p . The multi-state voter model has been considered in the past 
theoretically in terms of mappings of coarsening of the Potts model on the Ising model 
with constant magnetization ( I15p or in terms of duality properties (16) and has found 
applications in understanding the fragmentation transition in adaptive networks ( FlTl) or 
in neutral models of biodiversity (fT8| IT9]) . Variant of the pure multi-state voter model, 
introducing non-equivalent states, have also been discussed in the literature ( ]2"Ut l2"Tt |2"2"]) . 




(1) 
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In this paper we focus in the study the ordering dynamics of the symmetric multi- 
state voter model, focusing in particular in the limit of a large number of initial states. 
Each agent can be in one of S different but dynamically equivalent states. Agents follow 
the same dynamical update rules than in the binary version, with time being update at 
every dynamical step as t — » t + 1/N. Generalizing the formalism of the binary S = 2 
case, we obtain exact mean-field expressions for the exit probability and the consensus 
time for the general case of S states. We find that the consensus time increases very 
slowly with the number of states, and its difference with the binary case saturates as 
S — > oo; we rationalize this finding by comparison with the case of the mutant invasion 
in the ordinary two-state voter model. We then investigate the dynamics of the multi- 
state voter model on low dimensional lattices. Lacking of specific analytical insights, we 
compare the phenomenology numerically observed with the one exhibited by the mean- 
field case, and point out similarities and differences, focusing on the effect of the number 
of initial states and their configuration. The paper is structured as follows. Sec. [2] is 
devoted to the analysis of the mean-field multi-state voter model. Sec. [3] reports on 
numerical experiments concerning the low-dimensional case. Finally, in Sec. [4] present 
our conclusions. 

2. Mean- field analysis 

The analysis of the standard voter model at the mean-field level, based on the Fokker- 
Plank equation (1231) can be easily extended to the multi-state case (j2lj) . Let ft 
denote a generic configuration of the system (not unique) with voters in state i, 
n = {ni,ri2, . . . ,ns}, with a normalization X]j n « = N. The probability of finding the 
system in the configuration n at time t, P{n, t), evolves in terms of master equation that 
is defined by the transition rates w(n' — >■ ft) from the state n' to the state n. At each 
time step only one voter changes its state, consequently we can write a new configuration 
n' of a transition n — > n' as n' = n i+ j_ = {ni, . . . , rij — 1, . . . , rii + 1, . . . , ns} , being j 
and % the state of the voter before and after the transition, respectively. The transition 
rates n — > n i+ j^ and n — > n,i_j + are given by 



where A = 1/N is the natural microscopic time step of the model, while the transitions 
rates n i+ j_ — > ft and nj_ J+ — > ft have the form 



It is now possible to derive the associated master equation. Under the diffusion 
approximation (23), valid for large N, we consider the frequencies of the states x« —rii/N 
and we rescale the time by a factor 1/N, so that one time step t measures N updates 
of the voters. A generic configuration is therefore denoted by x = {xi,x 2 , ■ ■ ■ ,%s}, an d 
lies in the standard simplex S s = {x e M. s \ X]f x i = !}• The set of the S vertices of 
the simplex, Bs = {<? G M 5 |e* = 6ij, i = 1, ... ,S} is the absorbing boundary of the 




1 Tlj Tin 



(2) 




1 m + 1 rij — 1 



1 rii — 1 rij + 1 
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dynamics. We note that the constraint £V Xj — 1 reduces the number of independent 
variables from S to S — 1, so we can choose x# to be dependent from the others. 
Expanding the master equation in terms of 1/N we finally obtain, up to order N~ 2 , the 
final Fokker-Plank equation in continuous time (124"|) 

1 5-1 2 5-1 

dtP(x, t) = -J2% M 1 - *i) p & *)] - E ^ *)] > ( 4 ) 

where ft = d/dxi. 

The Fokker-Plank equation for the multi-state voter model at mean-field level does 
not have a drift term. This implies that the ensemble average density of each state (xj) is 
constant in time. This observation allows to extend the calculation of the exit probability 
in the standard voter model to general boundary conditions in the multi-state case. Let 
us define the generalized exit probability E A (x) as the probability that the system, 
starting in some random initial configuration x, orders in some configuration e 1 G A, 
being A an arbitrary subset of the absorbing boundary Bs- During the evolution of the 
system, the average densities are conserved. Let us define the quantity 4>a = J2i\e'eA Xi ^ 
which is also conserved. In the final consensus state, <p A will have a value 1 with 
probability E A (x), and a value with probability 1 — Ea(x). We hence obtain the 
generalized exit probability 

E A (x) = <f> A . (5) 

The consensus time for a given initial condition x is given, on the other hand, by 
the general equation (123|1 

S-l 5-1 

- N = - Xi)%T N {2) - 2 XiXjdidjTNix), (6) 

i=l j<i 

subject to the boundary conditions 

T N (xeB s ) = 0. (7) 

We can solve this equation by noting that the consensus time T^{x) has to be symmetric 
under any exchange Xj -H- Xj with i, j = 1, . . . S . Thus, we can impose the ansatz form 

5 

T N (x) = ^ f J r (x i ), (8) 

i=i 

where the function J-(x) is independent of S. Introducing this ansatz into Eq. ^ we 
obtain 

s s 

1 = -- Y, - Xi)%?(xi) = X: Q(x t ), (9) 
1=1 1=1 

where we have defined 

g(x) = -—x(l-x)F , (x). (10) 

From Eq. ^ and the normalization condition for Xj, se can see that the only possible 
values of G(x) are Q(x) = const = 1/S or Q{x) = x. Considering now the solution 
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for the S — 2 case in Eq. ([!]), we can see that the correct solution is given by the 
second case, which leads, after integration of Eq. (10), applying the boundary conditions 
jr(i) = jr(o) = 0, leads to 

s 

T N (x) = -Nj2( 1 -Xi)Ml-x i ). (11) 

i=l 

This solution generalizes the "entropic" form corresponding to the standard voter model, 
Eq. ([I]), recovering in a simpler way the formal result obtained in Ref . (1251) . 

From Eq. (11), we can analyze the behavior of the system in the limit of a large 
number of initial states. In particular, considering the homogeneous initial conditions 
Xj = 1/5, we have 

T%(S) = N(S-l)\n(-^-). (12) 



That is, as we could naively expect, the consensus time increases with the number of 
states allowed (the system is initially more disordered the larger the number of states 
and therefore requires more time to become ordered), but its growth is very slow and 
saturates in the limit S 3> 1. In fact, in the worst case scenario in a finite system, in 
which S = N, we have T§ (S = N) = N(N - 1) \n[N/(N -1)]->N in the limit of large 
N, being only a factor 1/ ln(2) ~ 1.44 larger that the binary case S = 2. This result can 
be rationalized considering that, when S = N, we are effectively describing an initial 
condition in which every agent has a different state, and the ordering arrives when one 
of this individuals manages to impose its state at the population level. It is therefore 
not surprising that we recover the N behavior observed in the binary model when the 
initial condition consists of a given state of one species in a population of individuals of 
the opposite state, and, crucially, only the runs in which the state of the mutant gets 
fixated are considered. From Eq. [12] we can obtain the form in which the saturation at 
large S is reached, namely, 

1 -5^o = 1 -^- 1 ) 1 »(^t)= ! ^ +0 ^ (13) 

where the last expression is asymptotically valid in the limit of large S. 

Another interesting property of the ordering dynamics of the multi-state voter 
model is the number of different states at time t, starting from an initial condition with 
S states. We define the number of surviving states as s(t) = J2i ^(^) / where Si(t) = 
if Xi(t) = 0, and Si(t) = 1 otherwise. To compute it at the mean-field level, we start 
by considering the average consensus time (Tjy(s)) for a random initial configuration 
x with s different states that can be computed by averaging the consensus time T N (x) 
over all the initial conditions in the simplex S s , 

(T N (s)) = jf dxT(x), (14) 

where |«S S | = i s the volume of the standard simplex S s . The integral in 

Eq. (14) can be computed using the variables a n = J2i x i> which respect the constraint 
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< o~i < o"2 < • • • < o"s-i < = 1) and noting that 

cr s _! log(a s _i)(icr s _i . . . / da x = / 5-1 log(a s _i)rfcr s _i = — . (15) 

Jo Jo v s ~ l Y- s v s ~ z >- 

From here it follows that 



(Tn(s)) = ^-Jl, (16) 
s 

Now, assuming that the average time to go from s + 1 to s states is AT = (T/v(s + 1)) — 
(T N (s)) — Ns~ 2 , for s > 1, we have 

l sm „ ^ + AT)-,(t) _ _ J ^ 

^ SW " AT "AT" AT [U) 

By integrating and inverting this relation we obtain that the number of surviving states 
s(t) starting with random initial conditions with s(0) = S states decays as 

s(t) = S (^] for t < N, (18) 



expression which is valid far from the ordering time of the system ~ N. In the case 
t ^> N, we expect s(t) ~ const in surviving runs; that is, averaging over dynamical 
realizations that have not reached consensus at time t. In this case, we will expect that 
s(t), averaged over all runs, will decay as the survival probability ( 126|) . Assuming the 
exponential form derived in Ref. (26) for the standard voter model, we will expect to 
observe 

s(t) ~ exp(-2t/A0 for t > N. (19) 

In Fig. [T] we check this prediction by means of numerical simulations of the multi-state 
voter model on a complete graph. The plot shows the behavior of s(t) for homogeneous 
initial conditions, which is fully compatible with the analytical predictions in Eqs. (18) 
and (pb. 



3. Numerical results in finite dimensional lattices 

In this section we present and discuss the results of numerical simulations of the multi- 
state voter model on lattice of dimension d — 1 and d = 2, comparing them with the 
analytical results obtained in at the mean-field level. 



3.1. Consensus time 

We focus in the first place on the behavior of the consensus time T^{x) as a function 
of the initial densities of the different states. We consider the simplest case S = 3, 
parametrizing the initial configuration as x = {xi,X2,xs} = {x, a(l — x), (1 — a)(l — x)}, 
with x G [0, 1] and a G [0, 1]. This parametrization preserves normalization, Yli x i = 1> 
and has the advantage that, for a given value of a, the whole range of values of x can 
be explored. In Fig. [2] we plot the consensus time T N (a,x) computed in lattices of 
dimension d = 1 and d = 2 as a function of x, and for different fixed values of a. In 
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Figure 1. Number of surviving states s(t) as a function of time in the multi-state 
voter model on a complete graph of size N = 400, starting with homogeneous condition 
and s(0) = S = N states. We compare with eq. ( [18] ). In the inset we show s(t) for 
t S> N, averaged over all runs, compared with eq. fl9| ), while s(t) ~ const if averaged 
only over surviving runs. 
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Figure 2. Normalized consensus time Tn(x)/Tn(0) as a function of x for the initial 
configuration x = {x, a(l — x), (1 — Qt)(l — x)} in regular lattices of dimension d = 1 
(left) and d — 2 (right) of size N — 400 sites, compared with the analytical mean-field 
prediction Eq. |TI| ). 



order to get rid of possible numerical pre-factors in the consensus time, we normalized 
it by its value at x — 0, which takes the form Tjy(at) = —N[a ln(a) + (1 — a) ln(l — a)} at 
mean-field level. The numerical simulations for d = 2 fit quite precisely the theoretical 
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Figure 3. Rescaled consensus time TS (S) /Tff (S = N) starting from homogenous 
initial conditions, as a function of the number of states S/N in dimension d — 1 and 
d = 2, for N — 10 3 , compared with the theoretical mean-field prediction Eq. (12 1. 
In the inset we plot the quantity 1 — T$ (S)/T§ (S — N), showing the power-law 
decay with S. Error bars represent the standard deviation error on the average of the 
distribution. Each point is averaged over 10 5 runs. 



mean-field prediction developed in Sec. [2j with only slight deviations for small a and 
close to x ~ 0.5. On the other hand, strong deviations are noticeable in dimension 
d = 1, specially for small values of a. This result is in agreement with the expectation 
for the standard voter model, in which the consensus time Eq. ([T]) is expected to be 
exact only for d > 2 (jSj). 



3.2. Effect of the number of states 

We have seen that, at the mean field level, and for homogeneous initial conditions 
Xi = 1/S for i = 1, . . . , S, the consensus time increases with 5 towards it limit value 



T$(S = N) as with a power-law form, as given by Eq. (13). In Fig. [3] we plot the 
rescaled consensus time T$ (S)/T$ (S = N) as a function of S for lattices of dimension 
d = 1 and d = 2 and fixed size N = 10 3 . From this figure we observe once again that the 
d = 2 behavior is well fitted by the mean-field prediction, while the d = 1 case shows 
deviations for small S. Interestingly, increasing the number of initial states S reduces 
the deviation from the mean-field theory, in a way that the behavior for S — > N is very 



well fitted by Eq. (13) 
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Figure 4. Surviving states s(t) as a function of rescaled time t/N Q ^ for ad=l and 
d = 2 lattice of N — 400 nodes, starting with homogeneous condition and s(0) = 5 = N 
states. For t <§; N e g, the number of surviving states decays as s(t) ~ t -1 / 2 for d = 1 
and s(t) ~ i _1 logi for d = 2 (inset). For large i, s(t) decays exponentially in both 
cases. 

3.3. Number of surviving states s(t) 

At the mean field level, the number of surviving states s(t), starting from maximally 
heterogeneous conditions S = N, decays as s(t) ~ St' 13 in the initial time regime, 
with an exponent (3 = 1, crossing over to a exponential decay at large times. In Fig. 
[2] we show the number of surviving states s(t) as a function of time corresponding to 
numerical simulations on d = 1 and d = 2 lattices. 

From the results of Fig. [4], it is clear that the decay at the tail of the density of 
surviving states is, as expected, exponential; in particular, we find that in the large time 
regime s(i) ~ exp(— 1.5t) for d = 1, while s(t) ~ exp(— 2t) for d = 2. For t <^n eff, 
on the other hand, the effect of dimensionality leads to a different, slower decay of s(t) 
than the mean-field prediction. In particular, we see that in d = 1, s(t) ~ Nt" 1 ^ 2 , while 
in d = 2, numerical data can be fitted to the form s(t) ~ Nt" 1 logt (shown in the inset 

of Fig. g. 

The origin of the slowing down in the decay of the number of surviving states in 
low dimensions can be attributed to the formation of spatial domains of sites in the 
same state, which have to annihilate diffusively in order to reach the consensus state. In 
this line, the behavior of the number of surviving states in d = 1 can be understood by 
means of a simple argument: At a given time t > 1, there will be a number of surviving 
states s(t). Assuming that sites with the same state form clusters, the activity will be 
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Figure 5. Consensus time T(S)/N e g as a function of the number of initial states S 
on a d = 1 lattice with N = 10 3 nodes with an initial configuration made of blocks 
of voters in the same state, of initial length N/S, compared with the consensus time 
obtained with random homogeneous initial conditions. In the inset we show that the 
quantity T C (S)/T H (S) — 1 goes to zero with a power law behavior with exponent -1. 
Error bars are obtained as in Figure [3j Each point is averaged over 10 5 runs. 



driven by the diffusive ffuctuation of the boundaries of those clusters, which will have 
a length £ oc t 1 ^ 2 . The number of different clusters will thus be s(t) oc N/£ ~ Nt^ 1 ^ 2 , 
recovering the observed time dependence. 



3.4- Effects of correlated initial configurations 

Considering the multi-state voter model on a finite lattice allows to check the effects 
of correlated initial configurations in the dynamical approach to the consensus state, 
effects which should be particularly important in one-dimensional lattices. We have 
thus simulated the multi-state voter model in a d = 1 lattice, stating from an initial 
configuration of S states arranged in S contiguous blocks on length N/S in a lattice 
of size N. In Fig. [5] we plot the consensus time for T^(S) in this correlated initial 
conditions as a function of S, comparing it with the consensus time starting with 
uncorrelated homogeneous initial conditions, T$(S). We find that the effect of starting 
with correlated initial condition strongly slows down the achievement of consensus. As 
expected, the difference between the consensus time with correlated and homogeneous 
initials conditions approaches zero for S — > iV with a behavior compatible with a power- 
law form of exponent —1, i.e. 

i^/^-iV 1 for S N. (20) 



T N (S) \N 
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4. Conclusions 

In this paper we have addressed the general scenario of the voter model in which the 
number of different states allowed in the model can be larger than two, and, in the 
thermodynamic limit, even unlimited. At the mean- field level, we have derived exact 
equations for the expression of the exit probability and consensus time, which generalizes 
naturally the 'entropic' form observed for the two-states case. We have highlighted that 
in the limit of S — > oo the ordering time is only 1/ In 2 times bigger than in the binary 
voter model, and with a simple analytic argument we have found the decay of the 
number of surviving states in times. Finally, we have studied numerically the behavior 
of the multi state voter model on 1— and 2— dimensional lattices, and compared the 
results with the binary case. The consensus time in the d = 2 case is well predicted by 
the mean-field theory, while the uni-dimensional case behaves differently. Remarkably, 
it increases with the number of initial states with a power-law form as predicted by 
the mean- field theory, for both d = 1 and d = 2. We have also addressed the effect 
of correlated initial conditions on the consensus time, finding out that the relevance of 
this effect decreases with the number of initial states with a power law behavior. In 
summary, our results show that the number of states is not a trivial parameter in the 
voter model, and it affects the overall dynamics in subtle ways. 
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